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Abstract 

CN 

C\ ■ We review the technique of heavy quark mass expansion of various operators 

made of heavy quark fields using a semiclassical approximation. It corresponds 
to an operator product expansion in the form of series in the inverse heavy quark 
mass. This technique applied recently to the axial current is used to estimate the 
charm content of the 77, rj/ mesons and the intrinsic charm contribution to the proton 
spin. The derivation of heavy quark mass expansion for (Q75Q} is given here in 
detail and the expansions of the scalar, vector and tensor current and of (QV^uQ) 
t^J- ■ (a contribution to the energy-momentum tensor) are presented as well. The obtained 

results are used to estimate the intrinsic charm contribution to various observables. 

(N 
O 
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^ 1 Introduction 

Nowdays it is established beyond any doubts that the naive picture of light hadrons as 
made of three constituent quarks (for baryons) or qq pairs of constituent quarks (for 
mesons) is not complete. The DIS experiments revealed the rich sea structure of the nu- 
cleon, these experiments showed in particular that a considerable portion of the nucleon 
spin is carried by the strange component of the nucleon sea. Furthemore there are exper- 
imental facts which seems to suggest that a non-vanishing nonperturbative component of 
intrinsic charm is present in light hadrons [p], [|. 

We address the problem of intrinsic charm content of light hadrons from the point of 
view of the heavy quark mass expansion. The cc pairs in light hadrons, due to parametri- 
cally large mass of charm quarks, can appear in a light hadron as virtual state whose life 
time is short, of order l/m c . The nonperturbative (with typical momenta below heavy 
quark mass m c ) gluon and light quark fluctuations are slowly varying from "point view" 
of virtual cc pair, hence the heavy quark mass expansion is equivalent to the semiclassical 
expansion. This expansion allows to rewrite operators made of heavy quarks in terms of 
light degrees of freedom (gluons and light quarks). For a detailed discussion of the heavy 
quark mass expansion see ref. 0. 

Let us note also that in absence of a direct probe of gluons the open charm production 
is considered as the main source of information about nucleon's gluon distributions. In 
hard leptoproduction heavy quarks are produced in the leading order via the photon-gluon 
fusion (PGF). The leading graph for PGF can be related directly to gluon distributions 
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if one assumes that there is no intrinsic charm content in the nucleon (no c(x),c(x) 
and no Ac(x),Ac(x) at normalzation point /x = m c ). However now there are many 
evidences that, in principle, there might be considerable intrinsic charm component in 
the nucleon wave function even at low normalization point. For reliable extraction of 
gluon distributions from open charm electroproduction experiments it is necessary to 
have quantitative estimates of the intrinsic charm content of the nucleon. 

This paper will be organized as follows: In the first part we present the calculation of 
the expectation value of heavy quark currents in the background of gluon fields using a 
semiclassical approximation. This corresponds to an expansion in the inverse of the heavy 
quark mass 

(Q^x)TQ(x)) = j:^ n) 



where T denotes the Lorentz structure of the current and the are local expressions 
of the field strength depending on T. In section [2.1| we review the large m expansion of 
the fermion determinant appearing in our definition of the expectation value. In section 
2~2| we then outline the expansion of color singlet currents in general before we present in 



section |2.2.1| the expansion of the axial current used in [|[] in full detail. The connection 
to the expectation value of the axial vector current using the axial anomaly equation and 
some general restriction coming from this equation are given in section 2.2.2| . As further 
examples we present the expansion of the scalar current in section |2.2.3| , the vector current 
(section |2.2.4j) and the tensor current (section |2.2.5|) , respectively. In section |2.2.6| we 
finally show the result of the expansion of {Q\x)V ^ y Q{x)) , appearing in the energy- 
momentum tensor of QCD. 

In the second part we discuss the calculation of intrinsic heavy quark content of light 
hadrons as an application of the heavy quark mass expansion. In the case of charm 
content of r]',r] mesons and intrinsic charm contributions to the proton spin we reduce 
the calculations of these quantities to matrix elements which are already known either 
phenomenologically or were computed previously. In other cases, like intrinsic charm 
contribution to the nucleon tensor charge and to energy momentum tensor, the problem 
is reduced to matrix elements of gluon operators which can be estimated using various 
nonperturbative methods in QCD: lattice calculation, QCD sum rule, theory of instanton 
vacuum. 



2 Heavy quark expansion of currents in the back- 
ground gluon and light quark fields 

The expectation value of a color-singlet quark current made of heavy quarks in the back- 
ground of gluon and light quark fields can be written after integration out heavy degrees 
of freedom as: 

(Q\x)TQ(x)) = detDtr C)7 (x|-r>) , (2) 



2 



Here T denotes an arbitrary Lorentz-structure. Note that all calculations will be per- 
formed in the euclidean space-time, so the QCD Dirac operator reads: 



D = i f + im (3) 
where the covariant derivative is defined as 

V, = (d, - (4) 

and m is the heavy quark mass. For the used conventions and the euclidization see the 
appendix. Eq.(|2]) can now be expanded in a power-series of the inverse heavy quark mass, 
1/m, under the assumption that the gradient of the background field strength is small 
compared to m. The expansion of determinant of the Dirac operator in eq.(|2]) has been 
calculated by a large number of authors, see e.g. f| §| 0]. We briefly review the calculation 
of the determinant following M since we use the result of this calculation as a check of 



the expansion of a scalar current of heavy quarks in Section 2.2.3. 



2.1 Expansion of the determinant 



The expansion of the determinant for heavy quarks in eq. (§) yields divergences of various 
types. Since most of these divergences are connected with the determinant of the free Dirac 
operator we normalize the determinant with that in zero external field. For the remaining 
infinity which can be related to the logarithmic renormalization of the coupling constant, 
we use the so-called ^-regularization. Using the identity 



detD = det£> f = (detZT^D) 5 
the normalized and regularized determinant can be written as follows: 



det£> 
det.Dn 



exp 



C-reg 



1, d M 2s 
lim 



2 s-o ds T(s) 



dtt s - L Tr 



(5) 



(6) 



where Dq denotes the Dirac operator in the absence of external gluon fields and M is 
the regulator mass. The functional trace denoted by Tr in eq.(|6|) can be calculated with 
respect to any complete set of states. For further calculations it is convinient to compute 
functional traces in the basis of plane waves, so that 



Tr 



-tD^D 



tr 



c,7 



tr, 



c,7 



d 4 x 
d 4 x 



d 4 k 

d 4 k 

N 1 



ikx 



-tD^D 



, ikx 



-tD*(d pl -+d IJ .+ ik ll )D{d„-+d u +ik v j 



1 . 



(7) 



The unity in eq.(^) points out that the operators here act on unity, so that ■ 1 = 0. 
The further calculations are straightforward: the expression in eq. (|7|) can be expanded in 
powers of the covariant derivative, integrated with respect to k and the Lorentz indices 
summed. 
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Since the explicit calculation is given in |7]] we present here only some useful formulae 
and the final result for the determinant up to order 0(l/m 2 ). The square of the Dirac 
operator in eq.(|7|) with all differentiation operators shifted, — > + ifc M , is given by 



D\d^ -^d„+ ik^)D(d u ^d u + ik u ) = -V 2 + ^-F - 2 ifcV + k 2 + 



m 



where we have used that 



pa 

-WW 



i [V M , V„f = d,Al - d v A% + f abc A b A c 



a 



(9) 
(10) 



with the notations aF = cr^^F^ and a^ v = + hv,7y] applied. Expanding the exponen- 
tial function in eq.(^) the functional trace then reads 



Tr 



Att 2 
+ t 



-tD^D 



tr. 



tr c . 



d x e 



-tm 



d 4 x 



t 2 



d A k 



-t(k 2 +m 2 ) 



(27T) 



E(-i; 



n=0 



TV. 



a 



F - 2iJfeV • 1 (11 



1. 



+ t u ( - V a V V a V - -V a V 2 V Q + -F a(3 F a/3 
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V 2 V 2 V 2 - — V Q V 2 V Q V 2 + V a V 2 V 2 V Q + V 2 V Q V 2 V 
36 



180 

+V Q V 2 V /3 V Q V /3 + V 2 V Q V^V a V^ + V a VpV 2 V p V a 

-rr (V a V /3 V a V 7 V /3 V 7 + V a V /3 V 7 V Q V /3 V 7 
yu 

+V Q V /3 V 7 V Q V 7 V /3 + V a V /3 V 7 V /3 V a V 7 + V a V^V 7 V/3V 7 V a ) 
+ g {FapFapV 2 + F a p\7 2 F a i3 + V 2 F a pF a p 

— V 7 F Q , /3 V 7 F Q/3 — V 7 F a «F Q/3 V 7 — FasV^FasV-y) 



--F aP F l(3 F ia ) + G(V 8 ) 



12) 



Here we used that all contributions with an odd number of fc's vanish whereas all other 
integrals with respect to k yield 



d A k 



rvii-, . . . rviif)^ C 



1 



(2tt) 4 " 4vr 2 
with S^...^ denoting all possible contractions: 



(2t) 



-(n+2) -tm 2 



Ml---A t 2n 5 



(13) 



exp 



■1 


a 2 - 








[2 a 




0M1 • • 




</>=0 



(14) 



After rearranging the terms into gauge invariants and taking also the part without external 
gluon fields in eq.@ into account, the determinant up to order 1/m 2 can be written as 
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follows: 



detD 
detD 



cxp 



C-reg 

i 1 

720^2^2 
11 1 



d 4 x 



1 , (M 2 \ „ „ 
In [ — | tr c F a pF a p 



48vr 2 



^ c F a pFp 1 F ia + : , — - — - tr, 



14407T 2 m 2 



tr, 



V 



7' 



14407T 2 m 2 

V a , F a/3 



V Q , -Fa/3 



V 7 , F 7/3 



3607T 2 m 



7/3 



384tt 2 m 2 
i 1 



d 2 ti c F^F lP ) 



exp 



48tt 2 



1 , /M 2 \ „ ^ 
In | — | tr c F a/3 F a p 



(15) 



70n 2 2 ^ T cF a pFp^F^ a 

72V7T Z m z 



1 1 

+ 120^ 2 ^ 2 ~ trc 



V 7 , F 7/3 



+ o 



rrr 



(16) 



Note that in the last step partial integration has been used with all total derivatives left 



out. The effective action S e g,E = — In detD which can be yielded from fllBD , rotated back 
to Minkowski space, corresponds exactly to the result of || [7fl: 



s t 



eflf,M 



where equation of motion terms, which vanish in pure Yang-Mills theory 
have been neglected. 



V a , F a p 



(17) 
= 



2.2 Expansion of heavy quark currents 

In order to expand tr^^l^rlx) in eq.(^|) in a series of the inverse heavy quark mass we 
can use eq.(Bf) to rewrite it as: 



tr c , 7 (x|— 1» 



tr 



c,7 



d A k 



- ikx 



tr, 



c,7 



d 4 ^ i 



DW 



D \ Te ikx 



(2tt) 4 k 2 + 



n=0 



A; 2 + m 2 



The expansion in eq.(^) is again justified for small gradients of the gluonic fields compared 
to the heavy quark mass m. Depending on the Lorentz structure V some of the integrals 
might be divergent and need to be regularized, we choose the dimensional regularization, 
since the integrals in eq.(|I~8"D can then be calculated using: 



d k k fll k fl2 . . . k^ 2ri 



{2n) d (k 2 + m 2 



1 T(n-m-d/2) / 1 \™- m ~<V 2 

o^ 1 ...^ m \—^j ■ (19) 



(Anf 2 



T (n) 2 m 
5 



The number of terms contributing to a given order of 1/m is reduced by the fact that 
terms containing an odd number of 7 matrices or an odd number of fc's vanish due to the 



trace over Lorentz-indices and the integration with respect to k. The expansion of eq. (IB") 
then is straightforward. The result of the expansion must be gauge invariant because 
we expand the gauge invariant operator. In order to obtain explicitely gauge invariant 
result for heavy quark mass expansion a number of helpful identities based on the Bianchi 
identity: 

V Q , F M ] + [V 7 , F af3 ] + [v^, F 7a ] = , (20) 

can be derived. They will be presented in the following sections. 

We want to illustrate some technical details of expanding heavy quark currents on the 
example of the pseudoscalar density and the divergency of the axial- vector current, which 
are related to each other by the axial anomaly. Another motivation to show detailed 
calculation for these cases is that recently confusing results for these cases were reported 
in the literature fl, pT0|1 . Further we present the result of the expansion of scalar, vector 
and tensor currents and of (QV^^Q) appearing in the energy-momentum tensor of QCD. 



2.2.1 The pseudoscalar density 

For T = 75 the expansion eq. (|18|) has the form: 



tr Ci0 ,(a;|— -y 5 |x) = -imtr 



c,7 




(2tt) 4 k 2 + m 



(21) 



Collecting all terms which contribute up to 0(1/ m?) one gets: 



tec,i{ x \jffo\ x ) 



imtr, 



c,7 



d 4 k 



(2tt) 4 L (F + m 2 ) 3 4 



1 1 „ T, 

-aFaFj 5 



1/1 1 1 1 

+ j -V 2 aFaF + -aFV 2 aF + -aFaFV 2 - -0F0F0F 1 7s 

(k 2 +w?) A V4 4 4 8 

(aFaFkVkV + aFkVaFkV + aFkVkVaF 



(k 2 + m 



2\ 5 



+WaFkVaF + WkVaFaF + kVaFaFkV)^ + O ( — ) (22) 



327r 2 m 



'■■ap-yS ^cFapFjg 



.„ o n&aB^tS cFpaFpB F™/j 



+ 



1927r 2 m 



FapF^V 2 — F a pV pF^V p + F Q/ gV 2 F 7 5 



VpF a/3 VpF 75 + V 2 F a/3 F 1 s — VpF a/3 F 7( 5Vp + O ( — - 



1 



77T 



(23) 



Here we have used eq.([L9|) for the integration over k and the following results for Dirac 
traces: 



Fa.pFyijX,T^ 



F a fjFjgtVj 



7a7/37 7 7575 



Faf}F^gF ei ptY^ 



4£ aPysFafiF^fi j 

^F a ^FygF ef pt r 7 

16 lEap^sFpaFppF^ 



7a7/3777a767^75 



(24) 
(25) 
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Using the following identities: 



V„,F 



pi "aft 



V p , F 7(5 



V„,F 



Vp, F 7 <5 



VpF a( gV p F 7( 5 + F a/3 V pF^V p — V pFapF^V p 
V 2 F a pFjs + F aj3 V 2 F l5 — 2V ' p F al3 V ' p F l5 , 
F a pV 2 F 7 s + Fa^FrysV 2 — 2F a pSJ p F^$V p 



eq. 



can be written as 



F a pV 2 F l5 , (26) 
(27) 
(28) 



tr C)7 (x|— 7s |x) 



1927r 2 m 3 



327r 2 m 

-a/375 tr c 



Vp, F aj g 



-F 7 5 + 



V„,F 



pi r Q.j3 



Vp, F 7 5 



+4iF oa F 0/3 F v<5 +0[^- 



From the Bianchi identity fl2"0|) we obtain the following relations^ 



pi 



Vp, F Q/ g 



-2i 



Fnr\ 1 F, 



pa j ^ p/3 



+ 



Vp,F p/3 



[V/3, 



V F 

v p> - 1 pa 



'a/375 tr c 



Vp,-Fp/3 



V a , F 7( j 



so that 



a/375 tr c 



Vp, F Q/ g 



-^75 : 

£ a/375 tr 



-2ie 



V, 



a/375 tr c 
Vp,F p/3 



FpaFppFyfi F n ftF nrv F* 



pf3"pa"-yS 



V/3, 



V F 

v pi - 1 - pa 



^75 



— 4 iFpaFppFryS + 2 V a , Vp, F p/ 3 



^a/3 7 5 (—4 i tr cFpaFppFys + 2<9 a tr c 



Vp,F p/ g 



-F75 
-F75 



On the other hands it yields 



'a/375 tr c 



Vp, 



Vp, F a/3 



F-yS + 



Vp, F Q/ g 



Vp, F 7<5 



'a/375 tr c 



Vp, F a/3 



(29) 

(30) 
(31) 



(32) 



75 



-ap^fsdp tr ( . 



Vp, F^ 



75 



^a/375^ptr c Vp, F a pFyS 

1 



-F 



a/3 



Vp, F 7 5 



-afl^sdp tr c 



Vp, F Q/ gF 7 5 



2 

1 2 

-£af3~{8d tr c F Q/ gF 7 5 . 



(33) 



So we finally end up with 

tr C)7 (x|— 75 |x) 



1 



327r 2 m 

l2 



a/375 tr c F Q( gF 7 5 



3847r 2 m 3 60/37,5,9 trc ^-^ 5+ 967r 2 m 3£:Q/375,9atrc 



Vp, F p/3 



^75 



O 



l6iT 2 m 



tr c F a/ 3F a/ 3 



1927r 2 m 3 



d 2 tr c F a sF a s + 



487r 2 m 3 



^atr c 



Vp,^p/3 



^a/3 + 



rrr 



(34) 



1 In the calculation of M the factor of 2 was missing in the first identity, which led to the wrong result. 



where we have introduced the common notation F a p = ^£ a p y 8F 7 s- 
2.2.2 The divergency of the axial vector current 

Instead of expanding the axial vector current j^{x) = Q\x)^ IJ rf^Q{x) in the way outlined, 
we can use that the divergence of the axial vector current is given by 

dJl = 2mQt 75 g _ ^ 2 F; V F; V , (35) 

where the first term contains the axial current and the second is the axial anomaly term 
which arises due to quantum effects. The expansion of the divergence of the axial vec- 
tor current in terms of the inverse of the heavy quark mass is therefore reduced to the 
expansion of the axial current, which we have already performed before. 



Further the axial anomaly equation (35) has some general properties, which can be 



used to check our result for the axial current: 

First the rhs of eq. (|35|) vanishes in the limit of infinite quark mass. This can be 
understood by the fact that the regulator mass cancels the physical mass in the infinite 
mass limit because of the different sign in the definition of the regulator. Therefore we 
expect the order 0(l/m) term in the expectation value of the axial current multiplied by 
2m exactly to cancel the anomaly term. Indeed equation (|34D gives: 

1 0(l/m) 



2mtr Ci7 (x|— 75 |x) = —ti c F^F^ v . (36) 



i 

D " J '~ / 87 2 

Second if one thinks of the expectation value of the axial vector current as a local 
large m expansion in the background of gluon fields 

WHE^ffO*) (37) 

then due to equation (|35|) the expectation value of the axial current in the large m ex- 
pansion is 

2mtr Ci7 (x|— 7s|x) = £ — d^xj%\x) . (38) 
u n m 

This means that terms appearing in the expansion of the axial current must be of the 
form of a total derivative. The order 0(l/m 3 ) term in equation ([34]) exactly obeys this 
form 

1 C(l/m 3 ) j 

D^ 5 ^ 967r 2 m 2 

R p = dp tr c F af3 F af3 + 4 tr c V a i Fquj F^ v . (40) 



2mtr C!7 (x|— 7 5 |x) = 2 d^R^ , (39) 



The term f abc F^ u F^ a F^ appearing falsely in the expansion of the axial current in || [J 
cannot be represented as a total derivative of a local expressionF] and therefore violates 
the general argument given above. 

2 A straightforward calculation for the instanton field shows that / d 4 a;/ o6c F° F% a F° ^ 0. But for 
dimensional reasons this nonvanishing contribution can be excluded from being generated by a surface 
term if the instanton field is taken in the regular gauge. Therefore the integrand cannot be a total 
derivative. 



The expectation value for the divergence of the axial vector current in the background 
of gluon fields finally reads up to order 0(l/m 4 ) 



96ir 2 m 

2.2.3 The scalar current 



-\- 4 tr c V Q , F au F^ u ) + o 



rrr 



(41) 



Following the steps outlined in the introduction to this section the expansion of a scalar 
current in series of the inverse heavy quark mass yields up to order (9 (1/m 3 ): 

d 4 k 1 



tr Cj7 (x| \x) 



tr, 



c,7 



(2 



7T 



— 1 



(4tt 



■ £2 _|_ m 2 

mf I 1 - - I dtr r l 



E 

?1=0 



'V 2 - |F + 2iW N 



k 2 



rrr 



i X 7 - J6 - im) ■ 1 



7i 



3607r 2 m 3 
i 



^cFapF^F^ — 2 — t^ 2 tr c F a/3 F Q/3 



28807r 2 m 



7207r 2 m 3 



11 tr, 



V/3, -F/3 7 



«7 



tr, 



V a , -Fa/3 



V 7 , F 7/3 



(42) 



The infinite constant term can be cancelled by substracting the expectation value of the 
scalar current for vanishing gluonic background fields. Our result eq. (|42|) coincides with 



that obtained in ref. [|ll|] if we neglect the total derivative terms which were ignored in 
ref. 0. 

Actually the result eq. ([4^) with the total derivative terms neglected (and hence that 
of ref. IjTTJ) can be easily obtained from the expansion of the determinant of the Dirac 
operator (|15D, since 



f 1 

J d 4 a;tr Cj7 (x| — \x) 



dm 
d 

dm 



iln(detD)) 
iTrlnD) . 



(43) 



Our expansion of the scalar current (fi2"|) is in agreement with the result for the determinat 
in equation fll5D. 



2.2.4 The vector current 

The heavy quark expansion of the vector current up to order 1/m 3 gives exactly zero 

^(x\^\x) = + o(J^ . (44) 

This result can be easily anticipated from the fact that the vector current is C— parity odd. 
This implies that the first operator contributing to heavy quark mass expansion should 
contain at least three gluon fields, additionaly the vector current conservation requires 
that this operator has the following structure VG 3 . From counting of dimensions we see 
that such operator can contribute only at 1/m 4 order. 
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2.2.5 The tensor current 

For the color singlet tensor current we find that the first non-vanishing order of the 
expansion is 0(l/m 3 ), yielding 

I i 1 

tr C)7 (x|-<v|:c) - — — 

FapFapFfiv + FavFpfiFaP — F afl Fp u F a/3 + O (^^j ■ (45) 



x tr, 



We note that the rhs of the above equation vanishes in the case of 577(2) gauge group. 
Actually one can show that the lhs of eq. ((45|) is identically zero in the case of SU(2) 
gauge group. Therefore the fact that rhs of eq. (|4"5| ) vanishes for SU(2) gauge group is a 
powerful check of our calculations. 

In order to prove that lhs of eq. ( |45|) is zero in the case of SU(2) gauge group we use 
the following transformation: 

G = Ct 2 , 

where C is charge conjugation matrix in Dirac spinor space and r 2 is color SU(2) matrix. 
Under this transformation we have: 

Gr a G- 1 = -r aT 
Ga^G 1 = — a^ u 
GDG- 1 = D T 

where T is the transposition operation. The lhs. of eq. ( f4"o] ) should be zero, since 

tr Cj7 (a;|-^(T^|x) = tVc^xlG^a^G'^x) 

T 

^c n {x\^a^ v \x) . (46) 



Nullification of the heavy quark mass expansion of the tensor current for 577(2) gauge 
group implies that lhs. of eq.fl45|) is zero if it is computed in the field of single instanton. 



2.2.6 Expansion of {QV^ V Q). 

The energy-momentum tensor of QCD can be written in Minkowski-space as 

T ,U = _ g ^ CQCD _ + ^ rf)^ ; (47) 

where (/zz/) denotes the symmetrization of the indices. The large m expansion of the (not 



symmetrized) last term in eq. fl47|) yields in Euclidean space 



. . 1 . . -2 i / 1 \~3 / d\ . 
tr c , 7 (x|-V, 7 ^) = J^{- 2 ) r (--j ^tr c l 
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r( — r ( 2 — - ] (——Sfu, tT c F a pF a p + - tr c F ai/ F afJi 



(4tt)2 \rn*/ \ i) 
+ „„„ n ^ c F a pFp,yFj a 



7i 



P — tr 



720vr 2 m 2 

Ii^^(lltr c [V a ,[V„F, 7 
1 1 / 

V/3, F aj g 



57607r 2 m 2 ^ "^«^<tf 



V 7 , F 7/ g 



2880tt 2 m 
-30itr 



V F 



+74 i tr c 
-26itr r 



V a , F a /3 



F QM + 14itr c 
F^ - 26 i tr c 



V/3,-^ 

V a , F Q/ g 



+22 i<9 2 tx c F Pv F PlJL - 3 i^Q, tr c F Q/3 F a/3 



-26itr c 
+40 i tr c 



V„,F, 



V^, F au V/3, F a /3 

V/3,F/3 M 



- 26itr, 
+ 4 i tr c 



Vj,, F au V/3, F a /3 



V w , F a p 



V M , F a /3 



( — 



and 



tr Ci7 (x|iv M 7^> 0(m °' m ~ 2) = tr c , 7 (x|V^7,|^) 0(m °' m " 2) 



The Lorentz trace of eq. (|48f) can be compared with the expansion of the scalar 
in eq. ( f4"2"| ) since 

^c,~f(A-^ vlA x ) = - mtT c, 7 (x\ — \x) . 



(48) 

(49) 
current 

(50) 



24vr 2 



For the trace of the rhs of eq. fl4"8|) we find 

tr Ci7 (x|lv,7,N)° (m( ' m " ) 1 
^cF a pF a ^Fp^ + 



7i 



3607r 2 m 2 

7r 2 rr) 2 \ 



720vr 2 m 2 



2880vr 2 m 2 

V/3, F/3 7 



^cF a pF a 
d 2 tr c F a/ 3F a/ 3 



V , F Q /3 



V 7 , F 7/ 3 



which exactly coincides with the rhs of eq. (^4) multiplied by (— m). Eq. (f49|) 
agrees with the expansion of the vector current beeing zero up to C(l/m 3 ) since 

tr^xlV^T^x) - tr Ci7 (a;|-^V /1 7 I/ |x) = <9 M Xx ca (x\^ v \x) . 



(51) 
further 

(52) 



3 Intrinsic heavy quarks in light hadrons 

In light hadron processes heavy quarks may give contributions only through virtual effects 
which are supressed by the mass of the heavy quarks. Especially for the charm quark, 
whose mass m c m 1.4 GeV is not too large, virtual processes nevertheless may give not 
negligible contributions. In this section we discuss the applications of heavy quark mass 
expansion obtained in the previous sections. 
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In the following sections we use the "perturbative" normalization for the gluon field 
strength G a = F*/ g and rotate all expressions to Minkowsky space, see the appendix. 



3.1 Intrinsic charm in 77 and vj 



For the decay of the 5-meson into 7/ and f^-mesons in [|] a mechanism with virtual charm 
quarks was suggested. In this approach the Cabbibo favored process b — > ccs is followed 
by the conversion of the cc pair directly into 77'. Its contribution to the decay amplitude 
is therefore direct depending on the "intrinsic charm" component of the rj- meson which 
is usually characterized by the matrix element 



(0|c7^7 5 c|r/(g)) = i/JV • 



(53) 



Using the heavy mass expansion of the divergence of the axial vector current ( pEl| ) , the 
constant f^' can be expressed up to the order l/m 2 c by 



Ac) 

Jri> 



1 a 
(01— G a G^'V) • 



(54) 



Here we neglected the term proportional to V a , G® in (41) which vanishes in pure Yang- 



Mills theory. We now can estimate the value of the constant f^, 



(c). 



In' ~ 



-2 MeV 



where we have used 



a 



— G a uy G llu,a \r]') = 0.056 GeV 3 



(55) 
(56) 



obtained in 12 



In QCD the omitted term can be related to the matrix element 



a s 

47T 



(57) 



q=u,d,i 



using the equation of motion. A rough order of magnitude estimate for the contribution 
of the omitted term (0) to f^, using the results of [[HJ indicates an deviation at the level 



of 0.3 MeV to the value (|55|). A more careful analysis of the omitted matrix element (|57j ) 
can be done by using the instanton methods developed in [T^| which already have been 
applied by to calculations of higher twist corrections to deep-inelastic scattering. 
Our estimated value for fjjP is consistent with the phenomenological analysis in [|I~5|j 

(c) 

where the authors dervived the bound —65 MeV < /V < 15 MeV from the analysis of 
777' transition form factors. From the analysis of (rj, 7]', ?7 c )-mixing in |l6j the small value 
= -(6.3 ± 0.6) MeV was derived, taking into account off-shellness effects in the cc 
component of rj' also, the value \fff\ ~ 2.4 MeV was found in 



I7| . Further our value 



for /j;, is in agreement with the phenomenological bound \f^'\ < 12 MeV, obtained in 



(c)| 



TE| , and corresponds to the result f^j, m —2.3 MeV presented in [B5|. In Ref. |T^[ the 
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divergence of the axial vector current was computed using the triangle graph for the axial 
anomaly with massive fermions, neglecting possible 1/m 2 contributions like 

cohered fib fic fKQ\ 

from higher order diagrams. Indeed our calculation shows that such "truly nonabelean" 
operators do not contribute to the order 1/m 2 and our result fl54l) therefore is exactly 
given by the first term of the expansion in 1/m 2 of the triangle graph 



The small value (p5| ) for fty implies that the b — > ccs mechanism does not play a 
major role in the B — > Krf decay mode. 

Bigger values of f!S due to the operator fl5"8|) in the expansion of the axial current 



up to order \jm\ have been given in f§, where ffi « (50 - 180) MeV and in f| with 
&) » -(12.3 - 18.4) MeV. These results have been used by a number of authors for 



the analysis of the charm content in noncharmed hadrons (see e.g ETI, ITDl YT\ E3]), but 



since the operator ( p8| ) violates general properties of the axial anomaly and it also does 
not appear in explicit calculations (see section |2.2.2| ), results relying on || [J should be 
reconsidered. 

Analogously we can immediately estimate the constant f( c ' characterizing the intrinsic 
charm contribution to the ^-meson. Using 

(0\^G a ^G^' a \7]} = 0.020 GeV 3 , (59) 



obtained in |E| we find 

4 c) ^-0.7MeV. (60) 

Since in the case of the rj meson the contribution of the omitted term (|57]) can be of the 
same order as ffl itself, the estimate (|60|) must be considered as a poor. 

3.2 Intrinsic charm contribution to the proton spin 

Another application of our result for the heavy quark mass expansion of the divergency 



of axial vector current fl4"I ) has been given in [pH[] . In this paper the authors have shown 
that the intrinsic charm contribution to the first moment of the spin structure function 
gi(x, Q 2 ) of the nucleon is small contrary to the result of EH], ID . 



In ref. it was proven that the forward matrix element of the axial current in the 
leading order of heavy quark mass expansion can be computed as: 

(N(p,X)\^ l5 c(0)\N(p,S)} = -^^(N(p,S)\R,(0)\N(p,S)}. (61) 

Here the current -R^O) is given by eq. ([40]). Note that the first term in does not con- 
tribute to the forward matrix element because of its gradient form, while the contribution 
of the second one is rewritten, by making use of the equation of motion, as matrix element 
of the operator 

ry _ ~ 

(N(p,S) \^ 5 c(0)\N(p,S)) = —J-{N(p,S)\g £ ^ nv G^ f \N(p,S)) 

c f=u,d,s 

" 2m%SjP, (62) 



127rm^ 
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(2) 

The parameter f s was determined before in calculations of the power corrections to the 



is 

,(2) 



first moment of the singlet part of g\. QCD-sum rule calculations gave 

f s f> = o.o9 m 

estimates using the renormalon approach led to fj?^ = ±0.02 [2^] and calculations in the 
instanton model of the QCD vacuum give a result very close to that of QCD sum rule 



Inserting these numbers we get finally for the charm axial constant of the nucleon the 
estimate 

^ = -^ 2)!! t--5-Kr 4 (63) 
yA 127T b m 2 c v 1 

with probably a 100 percent uncertainty. Note that this contribution is of non-perturbative 
origin (therefore we call it intrinsic), so that it is sensitive to large distances, as soon as 
the factorization scale is of order m r . 



3.3 Intrinsic charm contribution to the nucleon tensor charge 



Using the results of section [2.2.5| we can estimate the intrinsic charm contribution to the 
tensor charge of the nucleon. The tensor charge of the nucleon is defined as: 



(iV(p,5)|ca^75c(0)|iV(p,5)) = 2ig£\p II S u -p u S l 



(64) 



Using the result of the section |2.2.5| and the identity 



— e 



/iua/3 



a 



a 3 



we obtain for the charm contribution to the nucleon tensor charge the following result: 



9t 



967r 2 m\m 2 N 



xe 



Px S p (N(p, S)\ti c gl\G aP G^G^ + 2 G au G^G af3 } \N(p, S)) + O 



ml 



(65) 



The matrix element in the rhs of above equation can be rouhgly estimated in the instanton 
vacuum using method of |f[4 |. As discussed in section |2.2.5| the gluonic operator in rhs of 
eq. Q6"5D is identically zero on one instanton. Therefore the first nontrivial contribution one 
can get from instanton-anti-instanton pair. If we compare the expression for the charm 
contribution to the nucleon charge with that for axial charge we see that the charm 
contribution to the tensor charge is suppressed by additional power of m^jm c and one 
power of instanton packing fraction (7r 2 p 4 /^ 4 ), however the tensor charge is enchanced by 
one power of a s (m c ) []. This allows us to make a rough estimate for the charm contribution 
of the tensor charge: 



0) 
9t 



m N 



{N c 



7T 2 p 4 



a Am,) m. 



(c) 

9 a 



10" 



(66) 



3 Let us note that the expansion parameter in the heavy quark mass expansion is g s G/m c , because 
the non-perturbative gluon field strength G ~ l/g s (c/. instanton field). Therefore <? s (m) accompanied 
by gluon field strength is not counted as suppression. 
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Factors of N c are written in a way to reproduce the large N c behaviour of the matrix 
element and to account for the fact that the operator in rhs of eq. ( 65|) is identically zero 
at N c = 2. 



3.4 Intrinsic charm contribution to the nucleon momentum 

The charm contribution to the nucleon momentum can be defined as: 

i 



M 2 (c) (// 2 ) = f 1 dx B x B 
Jo 



c{x B ) + c(xb) 



2(P-n) 



■(N(P)\ci(n-V)c(0)\N(P)} , (67) 



where c(xb) is the charm parton distribution normalized at the scale fi, which is assumed 
to be fi ~ m c . The light cone vector n is arbitrary non-collinear to nucleon momentum 
P. 

Now we can use the result of eq. ( f48D in order to estimate the charm contribution to 
the nucleon momentum carried by intrinsic charm quarks. 



M 2 (c) (/i) 



2(P-n) 2 



' nJl ' ] ' 1 (N(P)\nWtr c G%G a »\N(P)) 



4tt (2-f 3 



+ i^ (jV(?)|nV tr c G au G^\N(P)} 



+0 f i ) • (68) 



mi 



In derivation of this expression we neglected terms which are proportional to [V a ,G a/ g] 
which are suppressed by one power of g s (/i) 2 . The first term in eq. (^) is divergent^ and 
actually is related to the mixing of quark and gluon operators. We can rewrite eq. (|68| ) 
as follows: 



+ 7^-r- 2 ,1^' \ N(P)\nV tv c G ar G^G">\N(P)) + O f -L ) . (69) 



Here the first term, which is proportional to the momentum fraction carried by gluons 
M\ accounts for extrinsic charm. Note that the coefficient in front of this term is 
exactly the leading anomalous dimension 7^ = 4/3 which accounts for mixing quark and 
gluon twist-2 operators under QCD evolution. The intrinsic charm contribution is given 
by the second term, so that we have estimates: 



M (c),int ri n Si c (//) = (N(P)\nW tr c G av G^\N(P)} + O f -L ) .(70) 



We see that the momentum fraction carried by intrinsic charm in the nucleon is related 
to the value of nucleon matrix element: 

(iV(P)Kn" ig s (fi) 3 ti c G au G^G al3 \N(P)) . 



We show only the most singular term 
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One can easily see that this matrix element in the theory of instanton vacuum [Il| is zero 
in one-instanton approximation, the same as matrix element 



(N(P)\nVg s (/i) 2 tT c G\G afl \N(P)) . 

Keeping in mind that for instanton field G ~ l/y g we can write: 

{N{P)\n»n v igsjriHTcGwGfrG^lNjP)) 
(N(P)\n»n»g s (riHr c G«G Q ^\N(P)} 



A 2 , (71) 



where A is parameter of the dimension of mass whose value can be obtained using various 
nonperturbative methods in QCD: lattice calculation, QCD sum rule, theory of instanton 
vacuum. Generically we expect that this mass parameter is of order of typical strong 
interaction scale A ~ 1 GeV. Now we can rewrite eq. ( [70|) in terms of this parameter and 
momentum fraction carried by gluons in the nucleon at scale fi ~ m c as: 

^ (C) ' mtrmS1C M = W^ Mf(/i) + ° (4) • (72) 

If we assume that A 2 =few GeV 2 than we get the estimate for the charm contribution to 
the nucleon momentum: 

M (c), intrinsic^ = few x 1Q -3 (73) 

We see that the heavy quark mass expansion of local currents allows us to reduce the 
problem of estimate of intrinsic charm content of the nucleon to the calculation of the 
ratio fl7ip . The latter ratio can be computed using various methods of nonperturbative 
QCD, probably the most promising would be a calculation of this ratio in lattice QCD. 

Recent analysis of refs p6| , |27fl gives for j\^ c )> mtrinsic values a t the level of fraction of 
percent what is in agrrement with our estimate (|73|) . 

Let us note that, since we performed the heavy quark mass expansion of heavy quark 
part of energy momentum tensor not neglecting total derivatives, one can compute also 
its non-forward nucleon matrix element. From the non- forward matrix element of energy 
momentum one can obtain the total angular momentum carried by intrinsic heavy quarks 



in the nucleon using Ji's sum rules fl28fl . The corresponding estimates we shall report 
elsewhere. 



4 Conclusions 

In this paper we have computed the heavy quark mass expansion of various local heavy 
quark currents. The details of the technique are illustrated on the example of heavy quark 
mass expansion of the pseudoscalar density Q^fsQ. This operator plays an important role 
in problems related to intrinsic charm contribution to the proton spin and to intrinsic 
charm content of 77, rf mesons. 

We corrected the mistakes in refs. || |J for heavy quark mass expansion of the operator 
• In these papers large intrinsic charm contribution to the proton spin and to 
intrinsic charm content of rj, rj' mesons was obtained due to contribution of the operator 
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cabcfia rib /"» 

presented in refs. 



f abc G a ^ v G b UQ G c a ^ which appeared in heavy quark mass expansion of the operator Q75Q 



||. We showed that coefficient in front of this operator is identically 
zero (the result which actually follows from general properties of the axial anomaly ||), 
so that the physical effects based on presence of the above operator discussed in refs. [§, 
TDL 



J3|] are absent. 

For the first time we presented the full results^ for heavy quark mass expansion of the 
operators QQ (to the order 1/m 3 ), (to the order 1/m 3 ), d^Q-y^Q (to the order 

1/m 2 ), Q'jfj.Q (to the order 1/m 3 ), Qa^Q (to the order 1/m 3 ), and Q^^S/uQ (to the order 
1/m 2 ). 

The results obtained for heavy quark mass expansion allowed us to estimate the intrin- 
sic charm content of 7/, 77 mesons as well the charm contribution to the proton spin, nucleon 
tensor charge and to the fraction of nucleon momentum carried by intrinsic charm. In 
the case of charm content of 7/, 77 mesons and intrinsic charm contributions to the proton 
spin we reduce the calculations of these quantities to matrix elements which are already 
known either phenomenologically or were computed previously. In other cases, like intrin- 
sic charm contribution to the nucleon tensor charge and to energy momentum tensor, the 
problem is reduced to matrix elements of gluon operators which can be estimated using 
various nonperturbative methods in QCD: lattice calculation, QCD sum rule, theory of 
instanton vacuum. 

We made here rough order of magnitude estimate of matrix elements of gluon operators 
appearing in heavy quark mass expansion of tensor current and of energy momentum 
tensor using instanton model of QCD vacuum. More quantitative estimates will be given 
elsewhere. The predictions for intrinsic charm contribution to various observables are 
summarized in Table 1. 



quantity 


our estimate 


Ac) 

Jr,' 


-2 MeV 


Ac) 


-0.7 MeV 


-"k 

9 A 


-5 • 10" 4 


(c) 

9t 


~ 10" 4 


J^( c ), intrinsic 


~ 10~ 3 



Table 1: Results for intrinsic charm contribution to various observables. 
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5 Appendix 

For the euclidization we use the follwing conventions: 

5 Not neglecting total derivatives and terms proportional to [V M , G p „] 
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l.r 



M ~~ " C 4, E j X M ~~ x k,E * d 4 XM 



i d 4 XE 



i^4,E , 



The covariant derivative therefore reads in Minkowski and in Euclidean space-time: 



M 



/x,E 



(0" - iA"(aO) M , 
(«9 M - L4 M (x)) E . 



The field strength, defined as 



F- v = d,A u - d v A, + r hc A\Al = i 



transforms as 



For the Dirac matrices we choose the conventions: 

7m = 74,e , 7m = i7fc,E , 
and 75 is defined within this paper as: 

7s,m = 7m = - i (7°7 1 7 2 7 3 ) M = (7i727s74) E = 75, 

With 



E • 



it yields 



.0123 
-M 



tr. 



~£()123,M — +1 — £l234,E 



757a7/37 7 75 



^a/3fS,E ■ 



The fermionic fields transform as 



^M = 4>E , lfo& = - l^E > 

so the Dirac operator in Euclidean space-time reads: 

D = i JI + im . 



(74) 



(75) 
(76) 



(77) 
(78) 

(79) 
(80) 



(821 



33) 



34) 



In section [2.1| we have used the following transformation properties for the effective action 
and the appearing operators: 



(FapF 



•Seff.M 

a /3> 



M 



(F a pF a p) E , {F a pF^ 1 F 



7 a 



M 



(FapFfrF, 



7"/E 



(85) 
(86) 
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